*Dear reader,

*If you have any questions about the dofile or the data, do not hesitate to contact me.

*Best regards,

*Rasmus Skytte

clear all
use "Data_from_experiments.dta"

***********************************************
**Creating indeces and other needed variables**
***********************************************

*creating dummy for party ID*
gen Dem0Rep1 = .
replace Dem0Rep1 = 0 if partyid == 5 | partyid == 6 | partyid == 7
replace Dem0Rep1 = 1 if partyid == 3 | partyid == 2 | partyid == 1
gen PID_undecided_equal_4 = partyid
replace PID_undecided_equal_4 = 4 if partyid == 8 | partyid == 9

*creating interaction terms manually*
gen civ_party = Dem0Rep1 * civility
gen pol_party = Dem0Rep1 * low_or_high
gen civ_pol = civility * low_or_high
gen civ_pol_party = civility * low_or_high * Dem0Rep1

*creating index for trust*
gen first_trust = (trust1 - 1)/3*(-1)+1 if trust1 != 5
gen second_trust = (trust2 - 1)/8*(-1)+1 if trust2 != 10
gen trust_index = (first_trust + second_trust)/2

*creating index for policy attitudes*
gen dril_attitude = (attitude_drilling - 1)/6*(-1)+1 if attitude_drilling != 8
gen air_attitude = (attitude_air - 1)/6*(-1)+1 if attitude_air != 8
gen att_index = (dril_attitude + air_attitude)/2

*creating index for affective polarization*
gen dematword = ((word_rating_democrats1) + (word_rating_democrats2*(-1)+6) + word_rating_democrats3 - 3)/4/3 if (word_rating_democrats1 != 6 & word_rating_democrats2 != 6 & word_rating_democrats3 != 6)
gen repatword = ((word_rating_republicans1) + (word_rating_republicans2*(-1)+6) + word_rating_republicans3 - 3)/4/3 if (word_rating_republicans1 != 6 & word_rating_republicans2 != 6 & word_rating_republicans3 != 6)
gen afpolword = repat - demat if Dem0 == 1
replace afpolword = demat - repat if Dem0 == 0
gen afpoltherm = (thermometer_rep - thermometer_dem)/100 if Dem0 == 1
replace afpoltherm = (thermometer_dem - thermometer_rep)/100 if Dem0 == 0
gen afpol_index = (afpoltherm + afpolword)/2
gen outpartyaffect = thermometer_dem/100/2 + dematword/2 if Dem0 == 1
replace outpartyaffect = thermometer_rep/100/2 + repatword/2 if Dem0 == 0
gen inpartyaffect = thermometer_rep/100/2 + repatword/2 if Dem0 == 1
replace inpartyaffect = thermometer_dem/100/2 + dematword/2 if Dem0 == 0

*creating variables for manipulation checks*
gen meanissuepol = ((perceived_issue_pol_drilling + perceived_issue_pol_air)/2-1)/3 if perceived_issue_pol_drilling != 5 & perceived_issue_pol_air != 5
replace meanissuepol = (perceived_issue_pol_drilling-1)/6 if perceived_issue_pol_drilling != 8 & sample == 3
gen meaninc = ((perceived_incivility_drilling + perceived_incivility_air)/2-1)/6*(-1)+1 if perceived_incivility_drilling != 8 & perceived_incivility_air != 8
replace meaninc = (perceived_incivility_drilling-1)/6*(-1)+1 if perceived_incivility_drilling != 8 & sample == 3

*labeling variables*

label variable low_or_high "High polarization (vs. low)"
label variable civility "Incivility (vs. civility)"
label variable Dem0Rep1 "Republican (vs. Democrat)"
label variable att_index "Policy attitudes (both issues)"
label variable afpol_index "Affective polarization (both issues)"
label variable trust_index "Trust (both issues)"
label variable inpartyaffect "Affect rating of in-party"
label variable outpartyaffect "Affect rating of out-party"
replace treatment_group = 5 if treatment_group == 0
label define treatment_group 1 "Low_pol_and_civ" 2 "Low_pol_and_inc" 3 "High_pol_and_civ" 4 "High_pol_and_inc" 5 "Control"
label values treatment_group treatment_group

*************************************************
**Making results and figure analysis in article**
*************************************************

***Scale reliabilities***

**Reliabilities reported in section called "Outcome measures"**
alpha first_trust second_trust if sample == 1 & passed == 1, std item
alpha afpoltherm afpolword if sample == 1 & passed == 1, std item

**Reliabilities reported in section called "Follow-up study 1: Using stronger incivility"**
alpha first_trust second_trust if sample == 2 & passed == 1, std item
alpha afpoltherm afpolword if sample == 2 & passed == 1, std item

**Reliabilities reported in section called "Follow-up study 1: Removing metacommentary"**
alpha afpoltherm afpolword if sample == 1 & passed == 1, item
alpha afpoltherm afpolword if sample == 2 & passed == 1, item

***Manipulation checks***

**Making Figure 2**
mean meaninc if passed == 1 & sample == 1, over(treatment)
estimates store R1
mean meanissuepol if passed == 1 & sample == 1, over(treatment)
estimates store R2

coefplot R1, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ci(off) yscale(range(0.2 0.8)) xlabel(,angle(45)) ////
ylabel(0.2 0.4 0.6 0.8) legend(off) yline(0.333 0.5 0.667, lpattern(dash)) yline(0.2 0.4 0.6 0.8, lpattern(dot)) ///
xlabel(1 "Low issue pol. / Civility" 2 "Low issue pol. / Incivility" 3 "High issue pol. / Civility" 4 "High issue pol. / Incivility") ///
ylabel(0.2 0.333 "Slightly polite" 0.4 0.5 "Neither rude nor polite" 0.6 0.667 "Slightly rude" 0.8, angle(0))

coefplot R2, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ci(off) yscale(range(0.6 1)) xlabel(,angle(45)) ////
ylabel(0.6 0.7 0.8 0.9 1.0) legend(off) yline(0.667 1, lpattern(dash)) yline(0.6 0.7 0.8 0.9, lpattern(dot)) ///
xlabel(1 "Low issue pol. / Civility" 2 "Low issue pol. / Incivility" 3 "High issue pol. / Civility" 4 "High issue pol. / Incivility") ///
ylabel(0.6 0.667 "Somewhat disagree" 0.7 0.8 0.9 1 "Strongly disagree", angle(0))

*Showing that debate was rated as "rude" in first follow-up as noted in section called "Follow-up study 1: Using stronger incivility"*
mean meaninc if passed == 1 & sample == 2, over(treatment)

***Main study***

**Making Figure 3 and results reported in section called "Results"**

*duplicating variables - only necessary to make the coefplot figures*
gen low_or_high2 = low_or_high_issue_pol
gen low_or_high3 = low_or_high_issue_pol
gen civility2 = civility_or_incivility
gen civility3 = civility_or_incivility
label variable low_or_high_i "High issue polarization (vs. low)"
label variable low_or_high2 "High issue polarization (vs. low)"
label variable civility2 "Incivility (vs. civility)"
label variable low_or_high3 "High issue polarization (vs. low)"
label variable civility3 "Incivility (vs. civility)"
label variable pol_party "High issue polarization (vs. low)"
label variable civ_party "Incivility (vs. civility)"

*results on trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if sample == 1 & passed == 1
estimates store Pooled
test low_or_high_issue_pol = civility_or_incivility

*results on attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party if sample == 1 & passed == 1
estimates store Pooled2

*results on affective polarization*
reg afpol_index low_or_high3 civility3 if sample == 1 & passed == 1
estimates store Pooled3

*Figure 3*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}") legend(off) ///
coeflabels(, labgap(2)) grid(none) format(%9.1g) mlabel mlabposition(12) xtitle("Change in scale points")

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy

gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy

gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color and thickness of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(width(medthick))) editcopy
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy

gr_edit plotregion1.plot4.style.editstyle area(linestyle(width(medthick))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy

gr_edit plotregion1.plot7.style.editstyle area(linestyle(width(medthick))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color and thickness of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(width(thick))) editcopy
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy

gr_edit plotregion1.plot5.style.editstyle area(linestyle(width(thick))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy

gr_edit plotregion1.plot8.style.editstyle area(linestyle(width(thick))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing size of value labels*
gr_edit plotregion1.plot3.style.editstyle label(textstyle(size(medium))) editcopy
gr_edit plotregion1.plot6.style.editstyle label(textstyle(size(medium))) editcopy
gr_edit plotregion1.plot9.style.editstyle label(textstyle(size(medium))) editcopy

*changing zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(pattern(dash)) editcopy
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

***Figure 4 and results for follow-up study 1 and 2 and 

*results on trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if sample == 2 & passed == 1
estimates store trust1
reg trust_index civility_or_incivility if sample == 3 & passed == 1
estimates store trust2

*results on attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party if sample == 2 & passed == 1
estimates store atpol1
reg dril_attitude Dem0Rep1 civility2 civ_party if sample == 3 & passed == 1
estimates store atpol2

*extra test - comparing the effect of attitude pol. in follow-up 2 on partisan attitude extremity*
gen part_att_extremity = dril_attitude if Dem0Rep1 == 1
replace part_att_extremity = (-1)*dril_attitude+1 if Dem0Rep1 == 0
reg part_att_extremity civility2 if sample == 3 & passed == 1

*results on affective polarization*
reg afpol_index low_or_high3 civility3 if sample == 2 & passed == 1
estimates store afpol1
reg afpol_index civility3 if sample == 3 & passed == 1
estimates store afpol2

*Figure 4*
coefplot (trust1 atpol1 afpol1), bylabel(First follow-up study) || (trust2 atpol2 afpol2), bylabel(Second follow-up study) byopts(xrescale) drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}") legend(off) ///
coeflabels(, notick labgap(2)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color and thickness of 90 percent intervals*
gr_edit plotregion1.plotregion1[2].plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plotregion1[2].plot1.style.editstyle area(linestyle(width(medthick))) editcopy

*changing color of markers*
gr_edit plotregion1.plotregion1[1].plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plotregion1[1].plot3.style.editstyle marker(fillcolor(black)) editcopy

*changing color and thickness of 95 percent intervals*
gr_edit plotregion1.plotregion1[1].plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plotregion1[1].plot2.style.editstyle area(linestyle(width(thick))) editcopy

*changing size of value labels*
gr_edit plotregion1.plotregion1[2].plot3.style.editstyle label(textstyle(size(medsmall))) editcopy

*changing zero line*
gr_edit plotregion1.plotregion1[1]._xylines[1].style.editstyle linestyle(color(black)) editcopy
gr_edit plotregion1.plotregion1[1]._xylines[1].style.editstyle linestyle(pattern(dash)) editcopy
gr_edit plotregion1.plotregion1[2]._xylines[1].style.editstyle linestyle(color(black)) editcopy
gr_edit plotregion1.plotregion1[2]._xylines[1].style.editstyle linestyle(pattern(dash)) editcopy

*changing headings*
gr_edit plotregion1.subtitle[1].style.editstyle fillcolor(none) editcopy
gr_edit plotregion1.subtitle[1].style.editstyle linestyle(color(none)) editcopy

*making labels for x-axis*
gr_edit plotregion1.xaxis1[1].edit_tick 6 0 `"Change in scale points"', custom tickset(major) editstyle(tickstyle(textgap(5)) )
gr_edit plotregion1.xaxis1[2].add_ticks 0 `"Change in scale points"', custom tickset(major) editstyle(tickstyle(textgap(5)) )

***Making Figur 5 (comparison with benchmark condition)***

gen T1 = trust_index
gen T2 = trust_index
gen T3 = trust_index
gen T4 = trust_index
gen T5 = trust_index
mean T1 if passed == 1 & treatment == 1
estimates store T1
mean T2 if passed == 1 & treatment == 2
estimates store T2
mean T3 if passed == 1 & treatment == 3
estimates store T3
mean T4 if passed == 1 & treatment == 4
estimates store T4
mean T5 if passed == 1 & treatment == 5
estimates store T5

gen De1 = Dem0
gen De2 = Dem0
gen De3 = Dem0
gen De4 = Dem0
gen De5 = Dem0
regress att_index De1 if treatment == 1 & passed == 1
estimates store R1
regress att_index De2 if treatment == 2 & passed == 1
estimates store R2
regress att_index De3 if treatment == 3 & passed == 1
estimates store R3
regress att_index De4 if treatment == 4 & passed == 1
estimates store R4
regress att_index De5 if treatment == 5 & passed == 1
estimates store R5

gen A1 = afpol_index
gen A2 = afpol_index
gen A3 = afpol_index
gen A4 = afpol_index
gen A5 = afpol_index
mean A1 if passed == 1 & treatment == 1
estimates store A1
mean A2 if passed == 1 & treatment == 2
estimates store A2
mean A3 if passed == 1 & treatment == 3
estimates store A3
mean A4 if passed == 1 & treatment == 4
estimates store A4
mean A5 if passed == 1 & treatment == 5
estimates store A5

coefplot T1 T2 T3 T4 T5 R1 R2 R3 R4 R5 A1 A2 A3 A4 A5, drop(_cons) vertical scheme(s1mono) fcolor(white) levels(95 90) yscale(range(0 0.50)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4 0.5) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(1) mlabcolor(black)) format(%9.2f) offset(0.0) lcolor(black) ////
xlabel(1 "Low pol. and civ." 2 "Low pol. and inc." 3 "High pol. and civ." 4 "High pol. and inc." 5 "Benchmark" ////
6 "Low pol. and civ." 7 "Low pol. and inc." 8 "High pol. and civ." 9 "High pol. and inc." 10 "Benchmark" ////
11 "Low pol. and civ." 12 "Low pol. and inc." 13 "High pol. and civ." 14 "High pol. and inc." 15 "Benchmark", labsize(small)) ////
xline(5.5 10.5) msymbol(O) mcolor(black) ciopt(color(black)) text(0.56 3 "Trust" 0.56 8 "Attitude polarization" 0.56 13 "Affective polarization")

gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot11.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot14.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot17.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot20.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot23.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot26.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot29.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot32.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot35.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot38.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot41.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot44.style.editstyle area(linestyle(color(black))) editcopy
gr_edit .style.editstyle margin(medlarge) editcopy

*********************************************
**Making results and figures for appendices**
*********************************************

***Appendix A5***

**Pooled results (excl. resp. who failed attention check)**

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1
estimates store Pooled

*attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if passed == 1
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

**Main study (excl. resp. who failed attention check)**

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & sample == 1 
estimates store Pooled

*attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1 & sample == 1 
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if passed == 1 & sample == 1 
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

**Follow up 1 (excl. resp. who failed attention check)**

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & sample == 2
estimates store Pooled

*attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1 & sample == 2
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if passed == 1 & sample == 2
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

**Follow up 2 (excl. resp. who failed attention check)**

*trust*
reg trust_index civility_or_incivility if passed == 1 & sample == 3
estimates store Pooled

*attitude polarization*
reg dril_at Dem0Rep1 civility2 civ_party if passed == 1 & sample == 3
estimates store Pooled2

*affective polarization*
reg afpol_index civility3 if passed == 1 & sample == 3
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(civility_or_incivility="{bf:Effects on trust in politicians}" civ_party="{bf:Effects on attitude polarization}" civility3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

***Appendix A6***

**Pooled results (incl. resp. who failed attention check)**

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility
estimates store Pooled

*attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

*Main study (incl. resp. who failed attention check)*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if sample == 1 
estimates store Pooled

*attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party if sample == 1 
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if sample == 1 
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

*Follow up 1 (incl. resp. who failed attention check)*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if sample == 2
estimates store Pooled

*attitude polarization*
reg att_index Dem0Rep1 civility2 low_or_high2 pol_party civ_party if sample == 2
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if sample == 2
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

*Follow up 2 (incl. resp. who failed attention check)*

*trust*
reg trust_index civility_or_incivility if sample == 3
estimates store Pooled

*attitude polarization*
reg dril_at Dem0Rep1 civility2 civ_party if sample == 3
estimates store Pooled2

*affective polarization*
reg afpol_index civility3 if sample == 3
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(civility_or_incivility="{bf:Effects on trust in politicians}" civ_party="{bf:Effects on attitude polarization}" civility3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) 

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

***Appendix A7***

**Note: Results from these analyses are also used in footnote 13**

*Pooled results*

reg inpartyaffect low_or_high_issue_pol civility_or_incivility if passed == 1 & (sample == 1 | sample == 2)
estimates store In_party_affect
reg outpartyaffect low_or_high_issue_pol civility_or_incivility if passed == 1 & (sample == 1 | sample == 2)
estimates store Out_party_affect

coefplot In_party_affect Out_party_affect, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Main study*

reg inpartyaffect low_or_high_issue_pol civility_or_incivility if passed == 1 & sample == 1
estimates store In_party_affect
reg outpartyaffect low_or_high_issue_pol civility_or_incivility if passed == 1 & sample == 1
estimates store Out_party_affect

coefplot In_party_affect Out_party_affect, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Follow-up study 1*

reg inpartyaffect low_or_high_issue_pol civility_or_incivility if passed == 1 & sample == 2
estimates store In_party_affect
reg outpartyaffect low_or_high_issue_pol civility_or_incivility if passed == 1 & sample == 2
estimates store Out_party_affect

coefplot In_party_affect Out_party_affect, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Follow-up study 2*

reg inpartyaffect civility_or_incivility if passed == 1 & sample == 3
estimates store In_party_affect
reg outpartyaffect civility_or_incivility if passed == 1 & sample == 3
estimates store Out_party_affect

coefplot In_party_affect Out_party_affect, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

***Appendix A8***

*Pooled results*

reg afpoltherm low_or_high3 civility3 if passed == 1 & (sample == 1 | sample == 2)
estimates store Thermomter_questions
reg afpolword low_or_high3 civility3 if passed == 1 & (sample == 1 | sample == 2)
estimates store Word_questions

coefplot Thermomter_questions Word_questions, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Main study*

reg afpoltherm low_or_high3 civility3 if passed == 1 & (sample == 1)
estimates store Thermomter_questions
reg afpolword low_or_high3 civility3 if passed == 1 & (sample == 1)
estimates store Word_questions

coefplot Thermomter_questions Word_questions, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Follow-up study 1*

reg afpoltherm low_or_high3 civility3 if passed == 1 & (sample == 2)
estimates store Thermomter_questions
reg afpolword low_or_high3 civility3 if passed == 1 & (sample == 2)
estimates store Word_questions

coefplot Thermomter_questions Word_questions, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Follow-up study 2*

reg afpoltherm civility3 if passed == 1 & (sample == 3)
estimates store Thermomter_questions
reg afpolword civility3 if passed == 1 & (sample == 3)
estimates store Word_questions

coefplot Thermomter_questions Word_questions, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

***Appendix A9***

*Pooled*
reg first_trust low_or_high3 civility3 if passed == 1 & (sample == 1 | sample == 2)
estimates store First_question
reg second_trust low_or_high3 civility3 if passed == 1 & (sample == 1 | sample == 2)
estimates store Second_question

coefplot First_question Second_question, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Main study*
reg first_trust low_or_high3 civility3 if passed == 1 & (sample == 1)
estimates store First_question
reg second_trust low_or_high3 civility3 if passed == 1 & (sample == 1)
estimates store Second_question

coefplot First_question Second_question, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Follow-up study 1*
reg first_trust low_or_high3 civility3 if passed == 1 & (sample == 2)
estimates store First_question
reg second_trust low_or_high3 civility3 if passed == 1 & (sample == 2)
estimates store Second_question

coefplot First_question Second_question, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

*Follow-up study 2*
reg first_trust civility3 if passed == 1 & (sample == 3)
estimates store First_question
reg second_trust civility3 if passed == 1 & (sample == 3)
estimates store Second_question

coefplot First_question Second_question, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) /// legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(col(1))

***Appendix A11***

**Intentions-competence**

*Intentions*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & (intentions == 1 & affective == 1 | intentions == 0 & affective == 0)
estimates store Pooled

*attitude polarization*
gen att_intentions = .
replace att_intentions = dril_attitude if intentions == 1
replace att_intentions = air_attitude if intentions == 0

reg att_intentions Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3  if passed == 1 & (intentions == 1 & affective == 0 | intentions == 0 & affective == 1)
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

*Competence*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & (intentions == 0 & affective == 1 | intentions == 1 & affective == 0)
estimates store Pooled

*attitude polarization*
gen att_competence = .
replace att_competence = dril_attitude if intentions == 0
replace att_competence = air_attitude if intentions == 1

reg att_competence Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3  if passed == 1 & (intentions == 1 & affective == 1 | intentions == 0 & affective == 0)
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

**Drilling-air traffic controllers**

*Drilling*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & affective == 0
estimates store Pooled

*attitude polarization*
reg dril_attitude Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if passed == 1 & affective == 1
estimates store Pooled3

*figure*

coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

*air traffic*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & affective == 1
estimates store Pooled

*attitude polarization*
reg air_attitude Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if passed == 1 & affective == 0
estimates store Pooled3

*figure*

coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

**Independents-partisans - what drives effects on trust?**

reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & partyid == 4
estimates store Pure_independents

reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & (partyid == 4 | partyid == 3 | partyid == 5)
estimates store Ind_incl_leaners

reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & (partyid == 1 | partyid == 2 | partyid == 3 | partyid == 5 | partyid == 6 | partyid == 7)
estimates store Partisans_incl_leaners

reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & (partyid == 1 | partyid == 2 | partyid == 6 | partyid == 7)
estimates store Partisans_excl_leaners

*figure*
coefplot Pure_independents Ind_incl_leaners Partisans_incl_leaners Partisans_excl_leaners, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12) legend(cols(1))

**Democrats vs. Republicans**

*Democrats*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & Dem0 == 0
estimates store Pooled

*attitude polarization*
reg att_index civility2 low_or_high2 if passed == 1 & Dem0 == 0
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if passed == 1 & Dem0 == 0
estimates store Pooled3

coefplot Pooled Pooled2 Pooled3, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" civility2="{bf:Effects on pol. attitudes}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

*Republicans*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & Dem0 == 1
estimates store Pooled

*attitude polarization*
reg att_index civility2 low_or_high2 if passed == 1 & Dem0 == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3 if passed == 1 & Dem0 == 1
estimates store Pooled3

coefplot Pooled Pooled2 Pooled3, drop(_cons) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" civility2="{bf:Effects on pol. attitudes}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

**DV after first or second issue**

*first*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & (order == 1 & affective == 1 | order == 0 & affective == 0)
estimates store Pooled

*attitude polarization*
gen att_first = .
replace att_first = dril_attitude if order == 1
replace att_first = air_attitude if order == 0

reg att_first Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3  if passed == 1 & (order == 1 & affective == 0 | order == 0 & affective == 1)
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

*second*

*trust*
reg trust_index low_or_high_issue_pol civility_or_incivility if passed == 1 & (order == 1 & affective == 0 | order == 0 & affective == 1)
estimates store Pooled

*attitude polarization*
gen att_second = .
replace att_second = dril_attitude if order == 0
replace att_second = air_attitude if order == 1

reg att_second Dem0Rep1 civility2 low_or_high2 pol_party civ_party if passed == 1
estimates store Pooled2

*affective polarization*
reg afpol_index low_or_high3 civility3  if passed == 1 & (order == 1 & affective == 1 | order == 0 & affective == 0)
estimates store Pooled3

*figure*
coefplot Pooled Pooled2 Pooled3, drop(_cons Dem0Rep1 civility2 low_or_high2) scheme(s1mono) xline(0, lcolor(gray)) msymbol(circle) ///
xscale(range(-0.23 0.23)) xlabel(-0.20 -0.10 0.00 0.10 0.20) levels(95 90) ///
headings(low_or_high_issue_pol="{bf:Effects on trust in politicians}" pol_party="{bf:Effects on attitude polarization}" low_or_high3="{bf:Effects on affective polarization}", labgap(-130)) legend(off) ///
yscale(alt) coeflabels(, notick labgap(-125)) graphregion(margin(l=65)) grid(none) format(%9.1g) mlabel mlabposition(12)

*changing color of markers*
gr_edit plotregion1.plot3.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot3.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot6.style.editstyle marker(linestyle(color(black))) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(fillcolor(black)) editcopy
gr_edit plotregion1.plot9.style.editstyle marker(linestyle(color(black))) editcopy

*changing color of 90 percent intervals*
gr_edit plotregion1.plot1.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot4.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot7.style.editstyle area(linestyle(color(black))) editcopy

*changing color of 95 percent intervals*
gr_edit plotregion1.plot2.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot5.style.editstyle area(linestyle(color(black))) editcopy
gr_edit plotregion1.plot8.style.editstyle area(linestyle(color(black))) editcopy

*changing color of zero line*
gr_edit plotregion1._xylines[1].style.editstyle linestyle(color(black)) editcopy

***Appendix A12***

*Trust*

mean trust_index if passed == 1, over(treatment_group)
estimates store R1

coefplot R1, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.40)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) lcolor(black)

mean trust_index if passed == 1 & sample == 1, over(treatment_group)
estimates store R1

coefplot R1, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.40)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) lcolor(black)

mean trust_index if passed == 1 & sample == 2, over(treatment_group)
estimates store R1

coefplot R1, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.40)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) lcolor(black)

*Attitude polarization*

regress att_index De1 if treatment == 1 & passed == 1
estimates store R1
regress att_index De2 if treatment == 2 & passed == 1
estimates store R2
regress att_index De3 if treatment == 3 & passed == 1
estimates store R3
regress att_index De4 if treatment == 4 & passed == 1
estimates store R4
regress att_index De5 if treatment == 5 & passed == 1
estimates store R5

coefplot R1 R2 R3 R4 R5, drop(_cons) vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.50)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4 0.5) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) offset(0.0) lcolor(black)

regress att_index De1 if treatment == 1 & passed == 1 & sample == 1
estimates store R1
regress att_index De2 if treatment == 2 & passed == 1 & sample == 1
estimates store R2
regress att_index De3 if treatment == 3 & passed == 1 & sample == 1
estimates store R3
regress att_index De4 if treatment == 4 & passed == 1 & sample == 1
estimates store R4
regress att_index De5 if treatment == 5 & passed == 1 & sample == 1
estimates store R5

regress att_index De5 if passed == 1& sample == 2

coefplot R1 R2 R3 R4 R5, drop(_cons) vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.50)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4 0.5) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) offset(0.0) lcolor(black)

regress att_index De1 if treatment == 1 & passed == 1 & sample == 2
estimates store R1
regress att_index De2 if treatment == 2 & passed == 1 & sample == 2
estimates store R2
regress att_index De3 if treatment == 3 & passed == 1 & sample == 2
estimates store R3
regress att_index De4 if treatment == 4 & passed == 1 & sample == 2
estimates store R4
regress att_index De5 if treatment == 5 & passed == 1 & sample == 2
estimates store R5

coefplot R1 R2 R3 R4 R5, drop(_cons) vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.50)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4 0.5) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) offset(0.0) lcolor(black)

*Affective polarization*

mean afpol_index if passed == 1, over(treatment_group)
estimates store R1

coefplot R1, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.40)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) lcolor(black)

mean afpol_index if passed == 1 & sample == 1, over(treatment_group)
estimates store R1

coefplot R1, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.50)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4 0.5) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) lcolor(black)

mean afpol_index if passed == 1 & sample == 2, over(treatment_group)
estimates store R1

coefplot R1, vertical recast(bar) barwidth(0.50) scheme(s1mono) fcolor(white) ciopts(recast(rcap) lcolor(black)) citop yscale(range(0 0.40)) xlabel(,angle(45)) ////
ylabel(0 0.1 0.2 0.3 0.4) legend(off) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) format(%9.2f) lcolor(black)